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The groundwave portion of atmospheric waveforms was examined to determine various 
characteristics of the radiation field from lightning discharges. Sixty-nine representative 
waveforms were selected from 1,029 atmospherics from 21 thunderstorms in the Oklahoma 
and North Texas area. The average amplitude and phase spectra, from 1 kc/s to 100 kc/s, 
are presented for several groups of atmospherics having distinguishable characteristics. 
Various relationships involving the total radiated energy, peak field strength, first half cycle 
length, spectral amplitude peak and frequency of spectra peak are presented. 

The "normal" type of atmospheric, composing 86 percent of the total, is predominately 
of positive initial polarity, has a spectrum peak near 5 kc/s and has well defined relationships 
between the parameters mentioned above. All "other" types of atmospherics are predomi- 
nately of negative initial polarity, have a spectrum peak from less than 1 kc/s to 18 kc/s and 
have no well defined relationships between the various parameters. 



1. Introduction 

The return strokes of a lightning discharge emit 
electromagnetic signals, called atmospherics, that 
usually attain a spectral peak within the VLF region, 
i.e., 3 to 30 kc/s, for the radiation field component. 
The spectra of these signals are very broad, extending 
from a few cycles per second to as high as 30 Mc/s. 
Aside from the interference or "noise" from lightning 
discharges, atmospherics have been utilized by many 
workers to study the mechanism of thunderstorm 
electricity, to detect and track thunderstorms, to 
study whistler mode propagation and to determine 
VLF and ELF propagation characteristics. It 
would be exhausting to properly credit all of these 
accomplishments and also extremely difficult because 
of much overlapping between the various categories. 
Examples of literature references covering these sub- 
jects are presented at the end of this paper. This 
paper will be limited to the presentation and descrip- 
tion of various characteristics of the radiation field 
from return stroke lightning discharges. The ground- 
wave portion of atmospherics from sources within a 
500 km range was used to minimize propagation 
effects. Corrections that were applied to the recorded 
data to overcome limitations of the recording equip- 
ment and alteration of the source function by ground- 
wave propagation characteristics are described. The 
data were separated into various groups by virtue of 
the wave shape and polarity of individual atmos- 
pherics. Characteristics of these groups are described 
in terms of the total radiated energy, peak field 
strength, first half cycle length, peak spectral ampli- 
tude, and frequency of the spectral peak. 



2. Collection of Data 

Atmospheric waveforms from a large number of 
lightning discharges in the Oklahoma and North 
Texas areas were recorded during the spring and 
summer of 1959. The results presented in this 
report are from the analyses of the groundwave 
portion of atmospherics recorded at a site near 
Lone Wolf, Okla. Selection of data was limited to 
observation periods occurring after sunset at this 
site to minimize interference from the first hop 
skywave. Other National Bureau of Standards 
recording stations were operated simultaneously 
with this station to determine the location of in- 
dividual lightning discharges for the present work 
and to record data for various VLF propagation 
studies. These other stations were located at 
Boulder, Colo.; Salt Lake City, Utah; Palo Alto, 
Calif.; and Maui, Hawaii. 

Each NBS station recorded the broadband vertical 
electric field of atmospherics, using a vertical an- 
tenna; direction of arrival, using a pair of vertical 
electrostatically shielded loop antennas with broad- 
band equipment; and the time of reception, using 
timing marks synchronized with the timing signals 
from stations WWV or WWVH. This information 
was displayed on multigun oscilloscopes and photo- 
graphed on continuously moving 35-mm film. The 
vertical and loop antenna channels included band- 
pass filters, wideband amplifiers and 24 jusec delay 
lines. Signals that exceeded a selected amplitude 
from the vertical antenna channel, prior to passing 
through the delay line were used to activate the 
oscilloscope displays. 
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The direction of arrival at each station was deter- 
mined from a conventional direction-finding (DF) 
display on one channel of a multigun oscilloscope. 
This DF channel was activated for approximately 
200 jusec at the beginning of each atmospheric. 
The direction of arrival was therefore determined 
from the initial portion of the signal, which is 
essentially vertically polarized, with an accuracy of 
about ±2° for each station. 

The time of reception was determined to an 
accuracy of about ± 1 msec. The limitation in tim- 
ing accuracy resulted from a relatively slow re- 
cording film speed of 3.3 in. /sec, and from the 
necessity of determining the correction for propa- 
gation time of the signals from the standard time 
stations to the recording sites. 

Atmospherics from individual lightning discharges 
were identified on the photographic records from 
each station by virtue of their times of reception 
and directions of arrival. The location of each 
discharge was determined by triangulation, using 
the direction of arrival indicated at each station. 

The recording sensitivities and triggering ampli- 
tudes were selected for each station to accommodate 
atmospherics originating within about a 500-km 
radius from the Lone Wolf site. The Lone Wolf 
equipment was adjusted to trigger from signals of 
plus and minus polarity with amplitudes equal to 
or greater than about 0.2v/m. The average dis- 
tance to the atmospheric sources selected for this 
study was approximately 350 km. 

The amplitude and phase as a function of frequency 
for the Lone Wolf vertical antenna channel are 
shown in figure 1. These curves represent the 
transfer characteristics for a transient signal 
presented on the cathode-ray tube of the multigun 
oscilloscope relative to the free space vertical 
electric field at the antenna. Calibration was 
performed b}^ inserting a voltage into a dummy 
antenna driving the cathode follower at the base of 
a 23-ft vertical antenna. The amplitude response 
was adjusted by high pass and low pass filters and 
sloped to 6 db points at 500 c/s and 100 kc/s. The 
phase response was zero from 4 kc/s to about 50 
kc/s, leading by 55 deg at 1 kc/s and leading by 
about 20 deg at 100 kc/s. This phase was computed 
for transient signals and therefore does not show 
the effect of an 8.6 degree phase lag per kilocycle 



1 1 1 1 1 1 — I r-i-j-n 

EFFECTIVE DELAY, 3.7 /aSEC— V 
MEASURED PHASE-^^^ 
PHASE RESPONSE ^--" 



EFFECTIVE PHASE- 




AMPLITUDE RESPONSE 



2 4 7 10 20 40 70 100 

FREQUENCY, kc/s 



Figure 1. Transfer characteristics of the vertical antenna 
channel. 



produced by the 24 /xsec delay line. The dashed 
curve shows the phase lag produced by the low pass 
filter and other electronic components. This approx- 
imates an effective delay of 3.7 /xsec shown by the 
dotted curve. It is the difference between the 
measured phase and the phase lag from a 3.7 /-isec 
effective delay that gives the effective phase response 
for a transient signal. 

3. Analyses 
3.1. Data Selection and Classification 

Twelve observation periods of 3 min each were 
selected for this study. The total number of 
atmospherics recorded from all distances was 2,763, 
or an average of one every 0.78 sec. Of this total, 
1,029 atmospherics, or one every 2.1 sec, originated 
from 25 thunderstorm areas within about 500 km 
of the Lone Wolf site. The average number of 
atmospherics per minute from each thunderstorm 
area was 4.6, with values ranging from 2.6 to 11.0. 
This latter value was obtained during severe activity 
centered some 300 to 400 km from the site. 

Examples of atmospheric waveforms are shown 
in figure 2. These are reproduced on the same time 
base but with individual amplitude scales. Dis- 
tances to the source are indicated and the first hop 
sky wave arrival time is shown by the arrow marked 
"s". Zero /xsec on the time base corresponds to the 
beginning of the VLF portion of the atmospheric; 
i.e., the return stroke of the discharge. The effect 
of the delay line used in the waveform signal channel 
of the recording equipment is evident in the initial 
part of each trace showing the background or 
noise preceding the main waveform. Some of these 
small deflections may be part of the radiation from 
a stepped leader or dart leader process. The 
higher frequency oscillations superimposed on the 
groundwave records are probably produced in part 
by variations in the return stroke current-moment. 
A number of workers have reported on these 
phenomena [Lutkin, 1939; Malan and Schonland, 
1951; Norinder, 1952; Pierce, 1955a and 1955b; 
Kitagawa and Brook, I960]. 

The groundwave and skywave portions of the 
waveforms shown in figure 2 are completely sepa- 
rated except for the waveform in figure 2b. The 
time interval between the reception of the ground- 
wave and first hop skywave pulses is a function of 
the distance to the source and the height of the 
ionospheric reflection layer, in this case, the D-iegion. 
The extent of interference between these pulses is 
determined by the time interval between reception 
and the pulse length of the groundwave. As a 
consequence of this, probably about 100 jusec of the 
second half cycle of the groundwave shown in figure 
2b was obscured. Prior to performing the spectral 
analyses of the groundwave pulses, under conditions 
represented by this sample waveform, the trailing 
edges were reconstructed as indicated by the dashed 
line. The differences that exist between the recon- 
structed and actual waveforms are believed to be 
small. Therefore, the error in the present analyses 
is not considered to be serious. 
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Figure 2. Examples of atmospheric waveforms. 



During the process of studying the waveform data 
of the recording site from sources within about 500 
km, it became evident that waveforms with per- 
sistent characteristics were observed during various 
observation periods. The author has selected six 
main groups into which all of the observed wave- 
forms could be placed. Each of these groups was 
subdivided into two groups determined by the 
initial polarity of each waveform. It should be 
emphasized that this grouping is not intended to be 
a new system of classification of atmospheric wave- 
forms. Some of these groups may be redundant 
with regard to the physical properties of lightning 
strokes. Atmospherics placed in groups 1 and 2 
might well be considered as a single "normal" type 
group. Also, the atmospherics placed in groups 3 
and 4 are similar in many respects, and could be 
considered as a single group. These groups were 
chosen only for the purposes of reporting the data 
and referencing within the paper. 

Group 1 -positive, represents those waveforms most 
often observed from lightning strokes. Representa- 
tives of this group are shown in figure 2a and 2b. 
This group is characterized by positive initial po- 
larity; an average zero axis crossing time or first half 
cycle length, of 55 jusec; a long second half cycle that 
is approximately 30 percent of the first half cycle 
amplitude; and often higher frequency perturbations 
along the waveform. Group 1 -negative, represented 
by the waveform in figure 2d, is characterized by 
negative initial polarity, and otherwise is similar to 
group 1 -positive. 

An example of group 2-positive is shown in figure 
2c. The second half cycle of waveforms in this 
group reaches a peak amplitude about GO percent 
that of the first half cycle. Otherwise, this group is 
similar to group 1-positive. Characteristics of 
waveforms in group 2-negative are the same as for 
group 2-positive except the initial polarity is negative. 

These first two main groups contain about 86 
percent of all atmospherics recorded at ranges of less 
than about 500 km. In the remaining analyses, 



atmospherics in these groups will be referred to as 
normal atmospherics. 

Characteristics of waveforms in groups 3, 4, and 5 
are similar in many respects. The samples shown in 
figure 2e, 2f, and 2g have first half cycle lengths 
much shorter than observed from normal atmos- 
pherics, and the waveforms are usually free of large 
higher frequency perturbations. Group 3-negative 
is characterized by an extremely small second half 
cycle amplitude. The second half cycle of group 
4-negative is about 50 percent that of the first half 
cycle. Group 5-negative is characterized In a very 
large second half cycle, about equal in amplitude to 
the first half cycle. The positive counterparts of 
groups 3, 4, and 5 are similar in shape except for a 
reversal in polarity. 

A summary of the number of atmospherics and 
the percentage of positive polarity waveforms for 
each group and for combinations of groups is pre- 
sented in table 1. 

Table 1. Number of atmospherics in each group. 

GROUP NO. TOTALS %+ TOTALS %+ TOTAL %+ 

1+ 78 



I- 



47 



831 



94 



2 + 


43 


2- 


9 


3 + 


5 


3- 


52 


4 + 


8 


4- 


66 


5 + 


3 


5- 


12 



131 



94 



-H029 



20 



Only 69 atmospheric waveforms were selected 
from the total of 1,029 waveforms recorded from 
atmospherics within a 500 km radius of Lone Wolf. 
1 1 was first of all necessary to limit the number to be 
used for spectral analysis because of the extensive 
labor required for the processing and scaling and 
the considerable expense for the electronic computer 
analysis. Therefore, the author decided to select 
not over 10 percent of the total for further analysis. 
The first requirement for each waveform was that it 
be free from obvious interference, such as signals 
from other lightning discharge processes and from 
the first hop skywave. Discharge processes such as 
stepped leaders, dart leaders, intercloud and intra- 
cloud discharges, etc., from relatively near (50 km) 
thunderstorms, produced a considerable amount of 
interference or background noise during some of the 
observations. First hop skywave interference was 
more severe at the greater ranges and therefore very 
few waveforms from lightning strokes at distances 
greater than about 400 km were selected. The elimi- 
nation of waveforms with a large amount of interfer- 
ence resulted in reducing the number of analyzable 
waveforms to about 18 percent of the total number. 
All remaining waveforms in groups 2,3,4, and 5 were 
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analyzed. To complete the selection, every fifth 
lightning discharge from each thunderstorm area 
were chosen. Then all waveforms with group 1 char- 
acteristics, radiated from the lightning strokes com- 
prising the chosen discharge, were analyzed. 

It is reasonable to conclude that the present selec- 
tion is representative of the 1,029 atmospherics re- 
corded for this study, and indeed representative of all 
atmospherics originating from this particular locale. 
Of course it must be realized that small changes in 
the relationships of the various parameters might 
result from the selection of waveforms other than 
those chosen for spectrum analysis in the present 
case. 

Another classification, group 6, is that of whistler 
producing atmospherics. These atmospherics are 
characterized by large amplitudes and a smooth os- 
cillatory nature after the first few hundred micro- 
seconds of the waveforms. [Helliwell, Jean, and Tay- 
lor, 1958; Hepburn, 1959; Norinder and Knudsen, 
1961.] Several waveforms of this type were recorded 
from sources within 500 km distance. These were 
not included in the present study because of the long- 
pulse length, low frequency characteristics of the 
source and interference from the first hop skywave. 

3.2. Spectral Computations 

The Fourier spectrum of the vertical electric field 
of an atmospheric can be represented by the integral: 



E(u,d)-- 



= f '<?(*,<*>•< 



l dt 



(1) 



where co is the angular frequency (2wf), d is the 
distance in kilometers, t is the time in seconds and 
the function G(t,d) is the instantaneous value of the 
field in volts per meter. Integration starts at time 
2 = and ends at t = r, since 

G(t, d)=0 for t<0 and for t> r. (2) 

Now (1) can be written in the form 

E^,d) = \E(^,d)\'e-^ ( ^ d) . (3) 

The function \E(w,d)\ is the amplitude spectrum and 
the function </>(co,d) is the phase spectrum as a func- 
tion of frequency (2irf) and distance. 

The complex spectral component of 69 atmos- 
pherics was calculated at intervals of 1 kc/s from 1 
to 30 kc/s, every 2 kc/s from 32 to 50 kc/s, and every 
5 kc/s from 55 to 100 kc/s, on an electronic computer. 
The spectrum of each atmospheric pulse was cor- 
rected for recording equipment response and normal- 
ized to a common distance of 1 km as a matter of 
convenience for other studies utilizing atmospherics. 

This normalization to the 1 km value not only 
included the inverse distance relationship for ampli- 
tude; i.e., \E(w,d)\ocl/d (this assumes a flat earth, 
but is sufficiently accurate at short ranges), but also 
included the effects of propagation on the amplitude 
and the phase of the signal. The amplitude and 
phase components of groundwave propagated signals 



have been computed by Johler, Kellar, and Walters 
[1956]; Wait and Howe [1956]. Applying the com- 
putations for amplitude and phase, reported in the 
above references, to the spectra of each atmospheric 
pulse, it was possible to approximate the true ampli- 
tude spectrum and phase spectrum of the radiation 
component at 1 km. The amplitude attenuation 
and phase variation for a given frequency is of course 
a function of distance and ground conductivity. It 
was assumed for this study that the effective ground 
conductivity in the Oklahoma and North Texas area 
was 10" 2 mhos/m [Jordan, 1950]. 

The "corrected" 1 km value of the amplitude 
normalization factor can be defined as the amplitude 
at 1 km divided by the amplitude at the propagation 
distance for each atmospheric and for each frequency 
of the amplitude spectrum. These amplitudes were 
derived from the amplitude versus distance curves of 
Wait and Howe [1956]. These theoretical curves 
include the effects of finite ground conductivity and 
of normal atmospheric refraction. The amplitude 
spectrum computed for each atmospheric was multi- 
plied by the corrected 1 km normalization factor to 
determine, at each frequency, the spectral ampli- 
tude of the radiation component at 1 km. 

An "uncorrected" 1 km value of the amplitude 
normalization factor can be defined simply as the 
propagation distance in kilometers. 

Differences between the amplitude spectra re- 
ferred to 1 km using corrected and uncorrected 
normalization factors were small for all frequencies 
(1 to 100 kc/s) at the shorter ranges (100 to 200 km), 
amounting to less than 15 percent of the uncorrected 
value. As the propagation distance increased, these 
differences became more pronounced. At a range of 
500 km, for example, the differences between the 
corrected and uncorrected values of the 100 kc/s com- 
ponent of the spectrum amounted to about 50 percent 
of the uncorrected amplitude spectrum value. 

As indicated in (1), the integration time started at 
time £ = 0. Time t was therefore the instant when 
energy first arrived at the recording site from the 
stroke. The propagation of energy corresponds to 
the group velocity. In free space, the group velocity 
and the phase velocity are equal to the universal 
constant, c; i.e., ^3X10 5 km/sec. The group and 
phase velocities, however, are less than c for ground- 
wave propagation. The group velocity, for purposes 
of this study, was assumed to be that of the 100 kc/s 
component of each waveform corresponding to the 
initial arrival of energy. Theoretical values for the 
group velocity, in the form of delay time (in ^sec) 
relative to free space, have been computed by Johler, 
Kellar, and Walters [1956]. Theoretical values for 
the phase velocity, in the form of degrees (phase lag) 
relative to free space, have been computed by Wait 
and Howe [1956]. 

The phase spectrum </>(w, d), for each atmospheric 
was corrected by subtracting the phase value result- 
ing from the phase velocity at each frequency less the 
group velocity at 100 kc/s. It should be noted that 
these phase values were a function of finite ground 
conductivity and distance as well as frequency. 
These corrections to the recorded phase spectra were 
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relatively small, amounting to less than 30 degrees 
for frequency components between 1 kc/s and 100 
kc/s and for distances from 100 km to 500 km. 

A relatively large error in the spectra of the re- 
corded waveforms resulted from the transfer charac- 
teristics of the equipment. Therefore, the spectra 
were corrected by the corresponding response pre- 
sented in figure 1. As an example: at 1 kc/s, the amp- 
litude was increased by 2 db and the phase increased 
by 55°; at 2 kc/s the amplitude was increased by 0.4 
db and the phase increased by 19°; and so on. 

The spectra presented in figures 3 through 7 are 
averages of various atmospheric groupings. The 
amplitude spectrum, |l?(co)|, in micro-volt-seconds 
per meter (/xvsec/m) versus frequency in kilocycles 
per second represents the amplitude of the radiation 
component of the electromagnetic field at a distance 
of 1 km. It should be noted that actual recordings 
at this distance would involve extremely large static 
field components. Also, it should be realized that 
radiated energy is proportional to the square of the 
spectral amplitude component at each frequency. 
The phase spectrum, <£(co), in degrees versus fre- 
quency represents the phase of the radiation com- 
ponent. Positive degrees indicate a phase lag. 
The electronic computer was coded such that the 
analysis of a sine wave, where t=0 corresponds 
to the time the sine wave crosses the zero axis 
moving in a positive direction, will result in a phase 
of zero degreees. Also, the analysis of a cosine 
wave, where £=0 corresponds to the time a sine 
wave reaches the maximum positive excursion, will 
result in a phase of minus 90 deg. 

Maximum deviation of the spectral values of 
individual atmospherics in each group is indicated 
by the dashed lines on either side of the average 
spectrum. These values were determined only 
at frequencies of 1, 2, 5, 10, 20, 50, and 100 kc/s. 
It should be noted that the maximum deviation of a 
particular group depends on the sample size or 
total number in the group. Therefore, it would be 
expected that the maximum deviation would in- 
crease somewhat with a larger sample size. 

Computing these values for the amplitude spectra 
involved normalizing the individual spectrum to 
the average total radiated energy value for each 
group. This was accomplished by multiplying the 
spectral components computed at the above fre- 
quencies by the quantity \E(u)\ i ™ l ^-\E(u))\* J ve . The 
\E(go)\ p is the value of the amplitude spectrum 
normalized to 1 km, at its maximum or peak value. 
The superscripts ind and ave denote respectively the 
spectral peak for individual atmospherics and for 
the average of the group as determined by the 
relation 



20,000 



E T =9X10-*\E(»)\\;' 



(4) 



where E T is total radiated energy. This relation, 
obtained from the data presented in figure 11, will 
be covered in greater detail later. The average 
total radiated energy for each group is indicated on 




2 4 7 10 20 40 70100 
FREQUENCY, kc/s 



each figure. 



Figure 3. Average amplitude spectra. 

(a) Group 1-positivo, avc ©r=61,000j. 

(b) Group l-m^;ilivi\ avc frY = 35,300 j. 

(c) Group 2-positive, avc A>=97,000 j. 



Maximum deviations as shown on the phase 
spectra did not involve any further computations 
or normalizations. Here again, the maximum devi- 
ations would tend to increase for a larger number of 
samples. 

The amplitude spectra in figure 3 are averages 
for the normal type atmospherics. The group 
1 -positive spectrum is shown in figure 3a, group 
1 -negative spectrum in figure 3b, and group 
2-positive in figure 3c. Representative waveforms of 
these were presented in figures 2a, 2d, and 2c, 
respectively. The frequency of the maximum ampli- 
tude of the group 1 -positive spectrum was at 5 kc/s. 
The spectra of all atmospherics in this group reached 
a peak between 4 kc/s and 7.5 kc/s. Similarly, the 
peak of the group 1 -negative spectrum was at 5.5 
kc/s, and the range of values within the group was 
between 5 kc/s and 6 kc/s; the peak of the group 
2-positive spectrum was at 5 kc/s, and the range of 
values within the group was between 4 kc/s and 
6 kc/s. Marked similarity was observed between 
the spectra of individual atmospherics in the normal 
groups and consequently between the group averages. 
Perturbations in the higher frequency components 
are probably produced by variation in the return 
stroke current mentioned earlier. 

The average amplitude spectra for group 3-minus, 
group 4-minus and group 5-minus are presented in 
figures 4a, 4b, and 4c and correspond respectively 
to the representative waveforms shown in figures 
2e, 2f , and 2g. The spectra for groups 3 and 4 are 
very similar and undoubtedly result from similar 
discharge characteristics. Although the frequency 
at which the average amplitude reaches a peak for 
these two groups is apparently less than 1 kc/s, the 
spectral peak of individual atmospherics extends to 
6.5 kc/s for group 3 and 5.5 kc/s for group 4. 

Presented in figure 4c, is the average amplitude 
spectrum for group 5-negative, represented by at- 
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Figure 4. Average amplitude spectra. 

(a) Group 3-negative, ave -Er=22,100 j. 

(b) Group 4-negative, ave Et= 6,600 j. 

(c) Group 5-negative, ave jEt=225,000 j. 
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Figure 5. Average phase spectra. 

(a) Group 1-positive. 

(b) Group 1-negative. 

(c) Group 2-positive. 



mospherics similar to the sample in figure 2g. 
Atmospherics in this group are usually of larger 
amplitude than other types. The average spectrum 
peaks at about 12 kc/s with peak spectral values of 
individual atmospherics in the group ranging from 
9 kc/s to 17.5 kc/s. The spectral amplitude, above 
20 kc/s, decreases approximately as the negative 1.5 
power of the frequency. This is greater than the 
rate of decrease observed for the other groups which 
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Figure 6. Average phase spectra. 

(a) Group 3-negative. 

(b) Group 4-negative. 
(e) Group 5-negative. 

is proportional to about the negative unity power 
of frequency. 

The phase spectra for the normal atmospherics 
are given in figure 5a for group 1-positive, 5b for 
group 1 -negative, and figure 5c for group 2-positive. 
Maximum deviations in phase of individual atmos- 
pherics are indicated by the dashed lines. The 
phase curves in figures 5a and 5c are for normal 
atmospherics with initial positive polarity. It 
should be noted that the phase is leading by ap- 
proximately 90° at the frequency for which the 
amplitude frequency reaches its maximum value, see 
figures 3a and 3c. Likewise, the phase in figure 5b, 
negative initial polarity, is lagging by approximately 
90° at the frequency of peak amplitude. The phase 
above 10 kc/s closely approximates a linear function 
of frequency. Computed from figures 5a, 5b, and 
5c, respectively, the phase lag per kilocycle is about 
1 deg, 3 deg, and 2 deg. At frequencies below 10 
kc/s the phase also approximates a linear function 
of frequency, with respective values of 10 deg, 10 
deg, and 15 deg per kilocycle. 

The phase spectra for groups 3,4, and 5, negative 
polarity atmospherics only, are presented in figures 
6a, 6b, and 6c. The corresponding amplitude 
spectra were given in figures 4a, 4b, and 4c. For 
groups 3 and 4, the phase is lagging by about 90° 
at 10 kc/s. Contrary to the phase and amplitude 
relationships of other groups, this does not correspond 
to the peak amplitude frequency. The phase above 
10 kc/s, for these two groups, closely approximates 
a logarithmic function of frequency (phase ^ 200 log 
of freq. in kc/s— 100). The phase below 10 kc/s 
tends to vary linearly with frequency, at a rate of 
about 12 deg per kilocycle. For group 5, see figure 
6c, the phase lags by about 90° at the peak amplitude 
frequency, as was the case for normal atmospherics 
of negative polarity. The variation of phase is 
approximately linear with frequency: about 7.5 deg 
per kilocycle at frequencies below 30 kc/s and about 
3.5 deg per kilocycle at frequencies above 30 kc/s. 
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Figuke 7. Average spectra, normal atmospherics, normalized 
to 1 km. 

(a) Amplitude | E{a>) \ . 
(i)) Phase 0(w). 

The average amplitude spectrum for the normal 
positive initial polarity atmospherics is presented 
in figure 7a. The amplitude component reaches a 
peak at 5 kc/s. At frequencies below the peak, 
the spectral amplitude is proportional to the 0.5 
power of the frequency; and at frequencies higher 
than the peak, the spectral amplitude decreases 
approximately as the negative unity power of 
frequency. 

Presented in figure 7b, is the average phase spec- 
trum for the normal positive polarity atmospherics. 
The phase at 5 kc/s, the peak spectral amplitude, is 
leading by about 90°. The variation of phase is 
approximately a linear function of frequency. At 
frequencies up to 10 kc/s, the phase increases, i.e., 
becomes more lagging, at a rate of about 12 deg per 
kilocycle; and above 10 kc/s the phase increases 
about 1.5 deg per kilocycle. 

The average amplitude spectrum of the normal 
negative polarity atmospherics should have the same 
characteristics as that for the positive polarity 
atmospherics. The average phase spectrum for the 
negative polarities should be also very similar to 
that for the positive polarities except that the phase 
will be increased by 180° at all frequencies. 

The average spectra presented in figure 7 were 
computed from the analyses of 47 selected normal 
type waveforms. They represent 86 percent of all 
atmospherics recorded at distances less than about 
500 km. 

3.3. Relationships 

There are a number of parameters involving the 
wave shape and the analyses of atmospherics that can 
be compared in an effort to better understand the 
phenomenon. It is the purpose of this section to 
present the relationships found to exist between 
various readily measurable parameters. 

The electromagnetic radiation from a lightning 
discharge contains a very small amount of the total 
available energy. Only about 1 part in 10"' or 10° 
is partitioned into radio energy. For the purpose al 
hand, it is desired to calculate the radio energy 
radiated into the earth's atmosphere by lightning 
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Figure 8. Pulse amplitude versus total radiated energy. 

discharges. The energy density, E d in joules per 
square meter, of the transient G(t, d) is represented 

by 

E *=iL£ G «> d)2dt (5) 

where G(t, d) and r were defined in (1) and (2). 
Assuming that the source can be represented by a 
vertical dipole, the radiation component of the 
vertical electric field is proportional to the cosine of 
the angle of elevation above the horizontal plane. 
The total energ}^ in the radiation term, E T , in joules, 
radiated into the atmosphere can be obtained by 
integrating the energy density, E d , over the surface 
of the radiation pattern of the source. Therefore, 



E 



-r c 

«/0 ./I) 



-AlZm- 2 E lt 



cos 3 6</d</<p 



(6) 



where r is the distance in meters from the source to 
the observation point, 6 is the angle of elevation 
measured from the horizontal, and <p is the angle 
measured in the horizontal plane. E T will be 
referred to as the total radiated energy in joules. 

There is a wide variation in the characteristics of 
individual groundwave pulses. Atmospherics from 
cloud-to-ground discharges usually consist of two 
half cycles. The first half cycle is larger in ampli- 
tude and shorter in duration than the corresponding 
values of the second half cycle. Variations in 
amplitude and duration of the first and second half 
cycles were the basis of selecting the various groups 
in section 3.1. In the following presentations, 
atmospherics in groups 1 and 2, the normal atmos- 
pherics, are represented by circles (O), atmospherics 
in groups 3 and 4 by dots (#), and group 5 atmos- 
pherics by crosses (X). The straight lines, drawn 
by eye through the data points, represent the rela- 
tionships between each of the two parameters for 
the normal atmospherics only. 

The distributions of the data points for the normal 
atmospherics in figures 8, 11, and 14 exhibit rela- 
tively good correlations. Fitting the line by least 
squares or some other more complex method would 
not necessarily increase usefulness of these data. 
The same is true for the lines drawn in figures 9 and 
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Figure 9. First half cycle time versus total radiated energy. 
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Figure 10. Pulse amplitude versus first half cycle time. 



10 although the data points are quite scattered. 
The accuracies of the lines presented in figures 12, 
13, and 15 are questionable because of the large 
scatter of data points. However, for purposes of 
this study, a line was drawn through the data points 
in figure 13, and then assuming this was correct, the 
relationships in figures 12 and 15 were computed. 
Thus, the relationship between the parameters in 
figure 12, for example, was computed from the 
relationships in figures 11 and 13. This was not 
necessarily the best method, but because of the large 
scatter of data in figures involving the spectral peak 
frequency, f , no substantial improvement in the 
significance or usefulness of these relationships would 
have been realized through the use of a more sophis- 
ticated method. A much larger sample of data 
would be necessary to improve the accuracies of 
these three relationships. Therefore, caution should 
be exercised in attempting to utilize the relation- 
ships involving f , presented here, for other related 
studies. 

The pulse amplitude versus total radiated energy 
is shown in figure 8. The pulse amplitude in volts 
per meter is the maximum instantaneous value of 
the radiation field, G(t) p , determined from each re- 




Figure 11. Peak spectral amplitude versus total radiated 
energy. 



corded waveform and normalized to 1 km. This 
normalization was performed by applying ground- 
wave propagation correction factors to the pulse 
amplitude depending on the propagation distance 
and the frequency peak of each spectrum. This is 
an approximation used for purposes of this report 
to relate the principal VLF components in the wave- 
form to a common distance. There is a close rela- 
tionship between these parameters as would be 
expected from (5) and (6) . The straight line drawn 
through the data points relates total radiated energy 
to the 2.2 power of the pulse amplitude or peak field 
strength of the radiation field at 1 km. Atmospherics 
with longer than average pulse lengths and with 
higher than average second half cycle amplitudes 
will plot to the right of the line drawn through the 
average of the normal type atmospherics. 

The relationship between the first half cycle time 
and the total radiated energy is shown in figure 9. 
The first half cycle time, t x in microseconds, often 
referred to as the length of the first half cycle, was 
determined by measuring the time between the ini- 
tial rise of the groundwave waveform and the point 
where the waveform crosses the zero axis. This 
point corresponds to approximately the time when 
the current moment reaches a maximum. The first 
half cycle time increases as the total radiated energy 
increases, but this is not a very close relationship. 
For other type atmospherics, represented by dots 
and crosses, the distribution on the graph is very 
irregular but in general the first half cycle time is 
much shorter for a corresponding total radiated 
energy than for the normal type. 

Figure 10 shows the relation between the pulse 
amplitude, G(t) p in volts per meter normalized to 
1 km, and the first half cycle time, t x in microseconds. 
The relationship is not very strong, but indicates 
the first half cycle time tends to increase as the 
pulse amplitude increases. The first half cycle 
times for other atmospherics are much less than for 
normal atmospherics at corresponding values of 
pulse amplitudes. 

The spectral peak versus total radiated energy is 
shown in figure 1 1 . This spectral peak is the maxi- 
mum value of the amplitude spectrum \E(o))\ p in 
microvolts seconds per meter, normalized to a 
distance of 1 km. The normalization includes 
groundwave attenuation as indicated in section 3.2. 
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Figure 12. Spectral peak frequency versus total radiated 
energy. 



A close relationship exists between these parameters 
as would be expected since 



.1: 



E T oc \E(a)\ 2 du. 



(7) 



The line through the normal type atmospherics re- 
lates total radiated energy to the 1.8th power of 
the maximum amplitude of the spectrum. Other 
type atmospherics tend to have greater total radi- 
ated energy than the normal type atmospherics for 
a corresponding spectral amplitude. This is under- 
standable since the other type atmospherics are 
characterized by a shorter pulse length, figure 9, 
and therefore a broader spectrum than the normal 
type atmospherics. A broad spectrum contains 
greater energy than a narrow spectrum for a cor- 
responding value of the maximum spectral amplitude. 

The relationship between the spectral peak fre- 
quency and the total radiated energy is shown in 
figure 12. The spectral peak frequency, / in kilo- 
cycles per second, is the frequency at which the 
amplitude spectrum reaches the maximum value. 
The relation is very poor but indicates in a general 
way that the spectral peak frequency decreases as 
the total radiated energy increases for normal type 
atmospherics. Group 3 and 4 atmospherics, dots, 
have a lower spectral peak frequency, and group 5 
atmospherics, crosses, have a much higher spectral 
peak frequency than normal type atmospherics for 
a corresponding total radiated energy value. 

Figure 13 shows the relation between the spectral 
peak and the spectral peak frequency. The relation- 
ship between these parameters is not very strong, 
but indicates that the peak frequency of normal 
atmospherics tends to decrease as the spectral 
amplitude increases. The spectral peak frequencies 
are lower for atmospherics in groups 3 and 4, and 
are higher for atmospherics in group 5 than for the 
normal atmospherics group. 
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Figure 13. Peak spectral amplitude versus spectral peak 
frequency. 
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A reasonably good relationship is found between 
the pulse amplitude and the spectral peak as shown 
in figure 14. This is expected since a larger wave- 
form amplitude produces correspondingly a larger 
spectral amplitude as indicated in (1). 

The first half cycle time versus the spectral peak 
frequency is shown in figure 15. The relationship 
is very poor. Group 3 and 4 atmospherics, dots, are 
scattered over the lower left of the graph; i.e., short 
first half cycles times and low spectral peak fre- 
quencies. Atmospherics in group 5, crosses, are 
characterized by spectral peak frequencies higher 
than for the normal type atmospherics. The group 5 
data points are reasonably well situated along an 
imaginary line parallel to the line drawn through the 
normal atmospherics which indicated that the 
spectral peak frequency is approximately inversely 
proportional to the square root of the first half 
cycle time. 

It would be extremely difficult to determine the 
relationships indicated by the straight line drawn 
through the normal atmospherics data points as 
shown in figure 15. This line was computed from 
the relationships found from the data presented in 
figures 9 and 12 and by solving for/ spectral peak 
frequency, in terms of t x , first half cycle time. 

4. Concluding Remarks 

An analyses summary giving the average values 
of each parameter for the various groups is presented 
in table 2. 





Table 2. 


Analyses summary. 






TOTAL 


^T 


Fo, 


IE(cu)l P 


tx 


G(t) P 


GROUP 


ANALYZED 


JOULES 


kc/s 


/y v sec/m 


/; sec 


v/m 


1 


36 


59,100 


5 


13,100 


53 


490 


2 


II 


97,000 


5 


20,500 


53 


460 


3 


8 


22,100 


<l 


6300 


20 


210 


4 


4 


6,600 


<l 


4,800 


23 


160 


5 


10 


225,900 


12 


18,700 


28 


920 


142 


47 


69,100 


5 


14,800 


53 


483 


3&4 


12 


16,900 


<l 


5,600 


21 


192 



The normal type atmospherics, represented by 
groups 1 and 2, have an average total radiated 
energy of about (39,000 j and a "spectral peak at 
about 5 kc/s. These values represent about 86 
percent of all atmospherics within a distance of 
500 km of the recording site, and are the primary 
source of interference at close ranges; i.e., less than 
2,000 km. However, at great ranges, the dominant 
waveguide mode of propagation is characterized by 
minimum attenuation near 15 kc/s [Chapman and 
Macario, 1956; Obayashi, Fujii, and Kidokoro, 1959; 
Taylor and Lange, 1959; Taylor, 1960a, 1960b; 
Wait, 1957a, 1957b, 1958b]. At 5 kc/s, the attenua- 



tion may be 4 to 12 db/1000 km greater than at 
15 kc/s, depending on path conditions. Therefore 
these atmospherics will be greatly attenuated over 
long distances. 

Atmospherics in groups 3 and 4 are characterized 
by small total radiated energies, about 17,000 j, 
and a spectral peak near 1 kc/s or less. These 
atmospherics will be severely attenuated at great 
ranges. 

The group 5 atmospherics have an average 
spectral peak near 12 kc/s and total radiated energy 
of about 226,000 j. These atmospherics will prop- 
agate to great distances with much less attenuation 
than experienced by the other groups. From these 
results, it can be concluded that the group 5 type 
atmospherics will become statistically more impor- 
tant at great ranges. If amplitude threshold 
triggering is employed for recording atmospherics 
from great ranges, all of type 5 will likely trigger the 
equipment, many of types 1 and 2 will be below the 
trigger level, and possibly all of types 3 and 4 will 
be below the threshold. 

The eight relationships determined from the data 
presented in figures 8 through 15 are listed below for 
convenience: 



E T =5XlO- 2 G(t) 2 p -' 
£V=4X10- 13 e- 7 
t x =U.lG(t)° p - 2 * 
E T =9X10-*\E(a) 
E T =5Xl0 18 Jo 18 

y =i4[^7(co)|-°- 1 

\E(u)\,=8.9G(t)l* 

/ =48C - 52 



(8) 

(9) 

(10) 

(11) 
(12) 

(13) 
(14) 
(15) 



The results presented here are particularly im- 
portant in the various areas where atmospherics 
produce a very serious interference problem. It 
should be possible, with this information, to improve 
methods that will discriminate against interference 
from the radiation field of lightning discharges. 

The use of atmospherics in VLF propagation 
studies has become important in recent years because 
of the application of this band to communication, 
timing standards, and navigation. Atmospherics 
are the only readily available signals at frequencies 
below about 15 kc/s, the lowest frequency commonly 
available from continuous wave communication 
emissions; or below about 10 kc/s, the lowest fre- 
quency normally used in the Radux-Omega naviga- 
tional system. Many workers, particularly in 
England, have reported VLF propagation character- 
istics from single station recordings for example 
[Chapman and Macario, 1956; Chapman and Pierce, 
1957a; Morrison, 1953]. These workers relied upon 
an assumed constancy of the radiation spectrum of 
atmospherics. In view of the results presented here, 
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that practice is hardly justified and unless extreme 
care is exercised in the selection of data for such work , 
erroneous conclusions may result. 

Previous work of Taylor and Jean ]1959] utilized 
atmospherics originating in the mountainous regions 
of Colorado and Wyoming. The cloud-to-ground 
strokes, for that study, occurred over ground of poor 
conductivity, i.e., <r<10~ 3 , and at elevations of 
6,000 ft above sea level or higher. The average 
total radiated energy was about 27,000 j. The 
present study involved analyzing a larger number 
of atmospherics than was previously possible and 
using the waveforms recorded from cloud-to-ground 
strokes occurring in the Texas and Oklahoma area. 
In this case, the strokes were over ground of good 
conductivity, i.e., <x>10~ 2 , and at elevations of 
about 1,000 ft above sea level. The average total 
radiated energy for the "normal" atmospherics was 
found to be about 69,000 j. It is indicated from 
these studies that topographical factors, such as 
ground conductivity and elevation, do control, at 
least in part, the characteristics of atmospherics. 
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